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ABSTRACT 

We have used semi-numerical simulations of reionization to study the behaviour of the power 
spectrum of the EoR 21 -cm signal in both real and redshift space. We have considered two 
models of reionization, one which has homogeneous recombination (HR) and the other in- 
corporating inhomogeneous recombination (IR). Considering the large scales first, we find 
that the predictions of these two models are similar. Both the real space H I power spec- 
trum P r (k) and the monopole moment of the redshift space H I power spectrum Pq(1c), 
fall sharply to a minima as the neutral fraction declines from Shi = 1 to 0.8 in the early 
stages of reionization. As reionization proceeds, P r and Pfi subsequently rise to a maxima 
at i ~ 0.4 — 0.5, and then declines in the later stages of reionization. In the early stages 
of reionization ( xu i > 0.8) the quadrupole moment of the H I power spectrum has a value 
consistent with P£/Pq = 50/49 predicted by the linear theory of redshift space distortion. 
This ratio falls abruptly at in i = 0.7, and is negative with P£ / Pq ~ —0.5 through the subse- 
quent stages of reionization. The predictions of the HR and IR models, we find, differ at small 
and intermediate length-scales. It is possible to qualitatively interpret the results of the simu- 
lations in terms of the fluctuations in the matter distribution and the fluctuations in the neutral 
fraction which have power spectra and dimensionless cross-correlation Paa(^), Pxx(k) and 
TZ = Pax/VPaaPxx respectively. We find 7Z = —1 at large scales through all stages of 
reionization. This provides a simple picture where we are able to qualitatively interpret the 
behaviour of both the real space and redshift space power spectra at large scales with varying 
xh i entirely in terms of a just two quantities, namely 5h i and the ratio P X x/ Paa- 
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1 INTRODUCTION 

The epoch when the neutral hydrogen (H I) in the inter-galactic 
medium (IGM) was reionized by the first luminous sources, is one 
of the least known periods in the history of our universe. Observa- 
tions of the CMBR (Spergel et al., 2003; Page et al., 2007; Komatsu 
et al., 201 1 ; Larson et al., 201 1) and absorption spectra of high red- 
shift quasars (Becker et al., 2001; Fan et al., 2003; White et al., 
2003; Fan et al., 2006; Willott et al., 2007; Goto et al., 2011) sug- 
gest that the epoch of reionization (EoR) probably extended over 
the redshift range 6 < z < 15 (Fan et al., 2006; Choudhury & 
Ferrara, 2006; Alvarez et al., 2006; Mitra et al., 2011). However 
these observations are limited in their ability to shed light on many 
important questions regarding EoR. What are the major sources of 
reionization? What are the typical sizes and the topology of the ion- 
ized regions at different stages? Observations of redshifted 21 -cm 
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radiation from neutral hydrogen hold the promise to answer some 
of these questions. The brightness temperature of the redshifted 21- 
cm radiation directly probes the H I distribution at the epoch where 
the radiation originated. It is thus possible to track the entire reion- 
ization history as it gradually proceeds with redshift. The presently 
functioning low frequency radio telescopes GMRT 1 (Swamp et al., 
1991) and LOFAR 2 , the upcoming MWA 3 and 21CMA 4 , and the 
future SKA 5 all cover the frequency range relevant for the EoR 
21 -cm signal, and this is one of the major goals for most of these 
telescopes. It is therefore very important to have a good picture of 
the expected signal in order to make forecasts for and correctly in- 
terpret the future observations of the redshifted 21 -cm radiation. 

There has been a considerable amount of work towards simu- 
lating the expected EoR 21 -cm signal. In particular, there have been 
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numerical simulations which use ray-tracing to follow the propa- 
gation of ionization fronts in the IGM (Gnedin, 2000; Ciardi et al., 
2001; Ricotti et al., 2002; Razoumov et al., 2002; Maselli et al., 
2003; Sokasian et al, 2003; Iliev et al., 2006; Mellema et al., 2006; 
McQuinn et al., 2007; Trac & Cen, 2007; Semelin et al., 2007; 
Shin et al., 2008; Iliev et al., 2008; Shapiro et al., 2008; Thomas et 
al., 2009; Baek et al., 2009). Such simulations are computationally 
extremely expensive, and it is difficult to simulate large volumes, 
and to re-run the simulations considering different values of the 
simulation parameters. Semi-numerical simulations which consider 
the average photon density in place of a detailed ray-tracing analy- 
sis provide a computationally less expensive technique to simulate 
the EoR 21-cm signal (Furlanetto, Zaldarriaga & Hernquist, 2004; 
Mesinger & Furlanetto, 2007; Geil & Wyithe, 2008; Lidz et al., 
2009; Choudhury, Haehnelt & Regan, 2009; Alvarez et al., 2009; 
Santos et al., 2010; Mesinger et al., 201 1; Zahn et al., 201 1). 

The fluctuations in the brightness temperature of the red- 
shifted 21-cm radiation essentially trace the H I distribution dur- 
ing EoR. The redshift space distortion caused by peculiar velocities 
also plays an important role in shaping the redshifted 21-cm signal 
(Bharadwaj et al., 2001; Bharadwaj & Ali, 2004). In fact, we ex- 
pect the peculiar velocities to introduce an anisotropy in the three 
dimensional power spectrum of the EoR 21-cm signal (Barkana 
& Loeb, 2005; Bharadwaj & Ali, 2005; Wang & Hu, 2006), very 
similar to the characteristic anisotropy present in the galaxy power 
spectrum (Kaiser, 1987). Barkana & Loeb (2005) have proposed 
that it may be possible to use this anisotropy to separate the effect 
of the peculiar velocities from the other astrophysical information 
present in the 21-cm power spectrum. 

Until recently, most simulations of the EoR 21-cm signal have 
not considered the effect of redshift space distortions. Some of the 
earlier work Mellema et al. (2006) and Thomas et al. (2009) have 
considered this effect while generating H I maps through their sim- 
ulations, but have not studied it's implication on the statistical prop- 
erties like the power spectrum of the brightness temperature fluc- 
tuations. Recently, Santos et al. (2010) and Mesinger et al. (2011) 
have included the effect of redshift space distortions in an approx- 
imate, perturbative fashion in their semi-numerical simulation and 
used this to study it's implications on the redshifted 21-cm bright- 
ness temperature power spectrum at different stages of the reioniza- 
tion. In a very recent work Mao et al. (2012) discuss the methodol- 
ogy to implement redshift space distortion in numerical simulations 
of reionization, and used this to study the 21-cm brightness temper- 
ature power spectrum during EoR. 

Most of the earlier semi-numerical simulations (e.g. Furlan- 
etto, Zaldarriaga & Hernquist 2004; Mesinger & Furlanetto 2007; 
Geil & Wyithe 2008; Lidz et al. 2009; Alvarez et al. 2009; Santos et 
al. 2010; Mesinger et al. 201 1; Zahn et al. 201 1) have assumed spa- 
tially homogeneous recombination which predicts strictly inside- 
out reionization where the most dense regions ionize first, the ion- 
ization subsequently propagating to lower densities. However, there 
are observations which indicate exactly opposite picture at the end 
of reionization, where the high density regions remain neutral (due 
to self-shielding) and the low density regions are highly ionized. 
Choudhury, Haehnelt & Regan (2009) have attempted to make their 
semi-numerical simulation consistent with these observations by 
incorporating the fact that recombination occurs faster in high den- 
sity regions. In these simulations reionization is inside-out only in 
the early stages. However, self-shielded, high density clumps re- 
main neutral in the later stages of reionization when inhomoge- 
neous recombination is taken into account. In this paper we follow 
Choudhury, Haehnelt & Regan (2009) to develop a semi-numerical 



code to simulate reionization, with the further improvement that 
we incorporate the effect of redshift space distortion due to pecu- 
liar velocities. We have used these simulations to study the effect 
of peculiar velocities on the EoR 21-cm signal, both with homoge- 
neous recombination and with inhomogeneous recombination. 

In this paper we have used semi-numerical simulations to de- 
termine the EoR 21-cm signal at different stages of reionization, 
and used the H I power spectrum to quantify the statistical proper- 
ties of this signal. We have calculated Pjg , (k) the H I power spec- 
trum in real space and its redshift space counterpart P^ , (k), and 
compared these two to asses the effect of peculiar velocities. The 
anisotropy of the 21-cm signal, quantified through various angu- 
lar multipoles of P^ : (k), is a very useful tool to study the ef- 
fect of redshift space distortion. In particular, we have studied the 
monopole and quadrupole moments of Pfj , (k) in order to identify 
the features characteristic of redshift space distortion at different 
stages of reionization. To our knowledge, this anisotropy has not 
been quantified using simulations in any of the earlier studies. Fi- 
nally, we attempt to interpret the results of our simulations, and 
compare these against the predictions of the simple, linear model 
proposed by Barkana & Loeb (2005). 

Unless mentioned otherwise, throughout this paper we present 
results for the cosmological parameters h = 0.705, Q m = 0.2726, 
fi A = 0.726, D. b h 2 = 0.02266. 

A brief summary of the paper follows. In Section 2. we present 
the semi-numerical technique that we have used to simulate the 
EoR 21-cm signal including the effect of peculiar velocities. Sec- 
tion 3. contains a brief discussion of the model prediction for the 
effect of redshift space distortion. These were used as reference val- 
ues in presenting amd interpreting the results from our simulations 
in Section 4. Finally, we discuss our results and conclude in Section 
5. 



2 SIMULATING REDSHIFT SPACE DISTORTION 
DURING REIONIZATION 

We have used a semi-numerical simulation to generates the 
H I ionization map during reionization. The simulation essentially 
starts from the dark matter distribution at a given redshift, and uses 
this to identify the sources of ionizing photons. These sources, 
along with the assumption that the H I traces the dark matter are 
used to construct a snapshot of the H I ionization distribution. Our 
simulation is based on the formalism proposed by Choudhury, 
Haehnelt & Regan (2009). This uses an excursion-set formalism 
as introduced by Furlanetto, Zaldarriaga & Hernquist (2004). The 
semi-numerical simulation provides us with the ionization field in 
the real space i.e. without the redshift space distortion. We briefly 
discuss the semi-numerical method that we have used in this work 
to simulate the brightness temperature fluctuations of the redshifted 
21-cm emission from EoR. 

We have used a Particle Mesh N-body code to generate the 
dark matter distribution. The spatial and mass resolution of the N- 
body simulation should be adequate to correctly resolve all the ion- 
izing sources that one is going to adopt in this semi-numerical sim- 
ulation. It is currently believed that the stars residing in the galax- 
ies are the major source of photons to reionize the universe (Yan 
& Windhorst, 2004; Stiavelli et al., 2004; Bouwens et al., 2005; 
Fan et al., 2006; Choudhury & Ferrara, 2006). The presently ac- 
cepted models suggest that dark matter halos having a mass M > 
10 9 ft -1 Mq host the early galaxies that contribute to reionization. 
We thus include all dark matter halos of mass M > 10 h~ l M® in 
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our semi-numerical simulation. Assuming that at least 10 dark mat- 
ter particles are required to constitute the smallest halo, the N-body 
simulation is required to have a mass resolution < 10 s h^ 1 Mq. 

We have generated the dark matter distribution at z = 8 using 
the Particle Mesh N-body code. The volume of the simulation is 
constrained by the 16 Gigabytes of memory available in our com- 
puter. We perform our simulation in a periodic box of size 85.12 
Mpc (comoving) with 1216 3 grid points and 608 3 particles, with a 
mass resolution M par t = 7.275 x 10 7 Mq. 

We identify halos within the simulation box using a standard 
Friend-of-Friend algorithm (Davis et al,, 1985), with a fixed linking 
length 0.2 (in units of mean inter particle distance) and minimum 
dark matter halo mass = 10M par t . 

The relation between the ionizing luminosity of a galaxy and 
its properties is not well known from the observations. In the semi- 
numerical formalism adopted here, we assume that the ionizing lu- 
minosity from a galaxy is proportional to the mass of it's halo. The 
number of ionizing photons contributed by a halo of mass M is 
given by 



M 

iV 7 (M) = N ion 

ran 



(1) 



where mn is the mass of a hydrogen atom and Ni OJl is a dimen- 
sionless constant. The value of Ni 0n is tuned so as to achieve the 
desired mass-averaged neutral fraction in the simulation. The ion- 
izing photon field is estimated on a grid which has a resolution 8 
times coarser than the N-body simulation. 

In our simulation we assume that the baryons follow the dark 
matter distribution and we also assume that each N-body particle 
has same hydrogen mass Mh- In the semi-numerical formalism, 
a region is said to be ionized if the average number of photons 
reaching there exceeds the average neutral hydrogen density at that 
point. Before applying the ionization condition we assume that the 
entire hydrogen contained in each particle is completely neutral. 
Using this assumption we calculate the H I density field in the same 
grid where the photon field has been generated. The photon and the 
H I density at each grid point are compared and a neutral fraction 
xh i is assigned to the grid point depending on the ionization con- 
ditions discussed below. 

In this work we consider two different models of reionization. 
In one model the recombination rate is assumed to be homogeneous 
(HR) and independent of density throughout the IGM. In the other 
model we have considered a density dependent inhomogeneous re- 
combination (IR). 

We next discuss how we implement the effect of peculiar ve- 
locities on the ionization maps generated. We first consider the dark 
matter particles to each of which we have assigned a total hydrogen 
mass Mh - The ionization map provides us with a neutral fraction 
xh i at each grid point of the simulation. For the ith particle in the 
simulation, we have interpolated the neutral fraction from its eight 
nearest neighbouring grid points to determine the neutral fraction 
x H i at the particle's position. We use this to calculate the particle 
H I mass M H i as 



M Hl =x Hl M H 



(2) 



This provides us with the H I distribution and the peculiar velocity 
associated with each H I element. We now consider a distant ob- 
server located along the x axis, and use the x component of the pe- 
culiar velocity to determine the particle positions in redshift space 



x + 



(3) 



aH(a) 

where a and H(a) are the scale factor and Hubble parameter re 



spectively. Finally we have interpolated the H I distribution from 
the particles to the grid, and used this to generate the EoR 21 -cm 
signal. 



3 MODELING REDSHIFT SPACE DISTORTION 
DURING REIONIZATION 

Coherent inflows into overdense regions and outflows from under- 
dense regions appear as enhancements in the matter density fluc- 
tuations observed in redshift space. This introduces an anisotropy 
(Kaiser, 1987) in P s (k,[i) the redshift space matter power spec- 
trum 



P 3 (fc,/i) = (l + /x 2 ) P r (fc) 



(4) 



where P r (k) is the real space power spectrum and p — k ■ h/k is 
the cosine of the angle between the wave vector k and the unit vec- 
tor n along the line of sight (LoS). Here we have assumed Q m = 1 
throughout, which is reasonable at the high redshifts of our inter- 
est. It is convenient (Hamilton, 1992, 1998; Cole et al., 1995) to 
decompose the anisotropy using Legendre polynomials Vi (jti) as 



(5) 



where Pi s (k) are the different angular multipoles of P s (k,p). 
Under the linear approximation (eq. 4), only the first three 
even moments have non-zero values P^(k)/P r (k) — 28/15, 
P£(k)/P r (k) = 40/21 and Pi(k)/P r (k) = 8/35 which are 
constant independent of k. 

Peculiar velocities have a similar effect on the brightness tem- 
perature fluctuations AT}, of the 21 -cm H I radiation from the high 
redshift universe (Bharadwaj et al., 2001; Bharadwaj & Ali, 2004). 
Expressing the brightness temperature fluctuations as ATi,(x, z) = 
T(z) r\H i (x, z) and considering the EoR where the spin tempera- 
ture is much higher than the CMBR temperature (Ts T 7 ) we 
have (Bharadwaj & Ali, 2005) 



T)H I (X, Z) = 

pH 



1 - 



(1 + z) dv\\ 



H(z) dr 



(6) 



where all the quantities in the r.h.s. refer to the position and epoch 
where the H I emission originated and T(z) — 4.0mK(l + 

2 ) 2 (%ot) (if) tJtTT- Here PR 1 refers t0 the H I density (which 
varies from position to position), pn is the mean hydrogen den- 
sity, r is the comoving distance from the observer and vu is the 
radial component of the peculiar velocity. We may express pn i 
using pH 1 1 pa = x h i (1 + 5) where xh i is the hydrogen neutral 
fraction and <5 is the matter overdensity which we have assumed 
to be the same as the hydrogen overdensity at the the large length- 
scales of our interest. Further, we may express the neutral fraction 
as£Hi = (xh i ) v where (xh i )v m the volume averaged 

neutral fraction and S x is the contrast in the xh i distribution. Note 
that the volume averaged neutral fraction (xh i)v and the mass 
averaged neutral fraction xh i refer to the average of xh i and 
xh i (1 + S) respectively, Assuming that S x , S, (dvu/dr) < lwe 
drop all quadratic and higher terms involving 5 X , S and (dvu/dr), 
to obtain {xh i)v = xh i whereby it is possible to express i)h i 
in Fourier space as 



f) H i (k) = x H i [A x + (1 + p 2 )A] 



(7) 



where t)h i , A and A x are the Fourier transform of j)h i , 5 and 8 X 
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Figure 1. Simulated H I maps in real and redshift space with different schemes of reionizations for various values of Sjj i- In this figure we use HR, HR- 
RS, IR and IR-RS to denote homogeneous recombination, homogeneous recombination with redshift space distortion, inhomogeneous recombination and 
inhomogeneous recombination with redshift space distortion respectively. The line of sight is along the x-axis and the coordinates are in Mpc. 



respectively. This gives the redshift space H I power spectrum (of 
t]h i ) to be (Barkana & Loeb, 2005) 

Pl I1 {k,n)=x 2 Hl [P xx {k) + 2(1 + ^ 2 )P Ax (k) + (1 +v) 2 Paa 

(8) 

where P xx and Paa are the power spectra of A x and A respec- 
tively, and Pax is the cross power spectrum between A and A x . 
We recover the real space H I power spectrum P r (k) if we set 
H = in eq. (8). Here, and in the subsequent discussion, we drop 
the subscript H I for brevity of the symbols. Throughout this pa- 
per we have used the model presented above [eq. (8)] to interpret 



the H I power spectrum. We use A(k) — \J P xx (k) / Paa (k) and 
n{k) = P Ax (k)/^/P xx (k)PAA{k), where TZ(k) is the cross- 
correlation coefficient between the fields A(fc) and A x (k), to ex- 
press the results of this model. We have the H I power spectrum in 
real space as 

P r = xh, (1 + 2ATZ + A 2 ) Paa . (9) 
In this model only the first three even angular moments of the red- 




Figure 2. The H I power spectrum in real space (top) and in redshift space (bottom) with the different reionization schemes indicated. For clarity the error bars 
are only shown for the case where xh i = 1. Note that the error bars increase with decreasing neutral fraction, and they are quite large for ijj j = 0.1. 



shift space power spectrum have non-zero values 

Po = x 2 Hl + \An + A 2 ^j Paa , (10) 

PS = %■ (J^ + ^aii\paa, (ii) 

Pl=x 2 Hl (j^jP A A. (12) 

Finally, we note that this model is based on the assumption 8 X <C 1, 
and terms of the order of S X S which appear in eq. (6) are ignored. 
While this is possibly a reasonable assumption in the early stages of 
reionization (xh i ~ 1), we may expect significant deviations from 
this model in the late stages of the reionization where xh i is small 
and we expect S x to exhibit large fluctuation of order unity. Lidz 
et al. (2007) and Mao et al. (2012) have considered models which 
incorporate the non-linear effects, but these are rather complicated 
and we have not attempted using these models here. 



4 RESULTS 

We have simulated the H I distribution, in both real and redshift 
space, for xh i values starting from 1.0 to 0.1 with an interval of 
Axh i = 0.1. For each value of xh i we have used 12 independent 
realizations of the simulation to estimate the mean and la error 
of the H I power spectrum. The simulations are all carried out at 
a fixed redshift z — 8. Ideally, one should adopt a model for the 
evolution of xh i with z and simulate each neutral fraction at the 
appropriate redshift. This, however, makes the results dependent 
on the model for the redshift evolution of xh i and also makes the 
computations rather cumbersome. We have sidestepped this issue 



by keeping z fixed at a value z — 8 which is in the redshift range 
that is considered favourable for observing the EoR 21 -cm signal. 
We interpret the simulations for different values of xh i as repre- 
senting different stages of reionization. 

Figure 1 shows the H I distributions at different stages of 
reionization in both real and redshift space for a single realiza- 
tion of the simulation. We see that the H I distributions in the HR 
and the IR models are nearly indistinguishable in the early stages 
of reionization (xh i > 0.7). In the early stages reionization pro- 
ceeds inside-out in both of these models with the high density re- 
gions ionizing first and the low density regions ionizing later. The 
H I maps of the two models become quite distinct at the later stages 
of reionization (xh i < 0.5). In the IR model we find many small 
isolated neutral H I clumps distributed within the regions which 
are completely ionized in the HR model. These small H I clumps 
correspond to the high density regions which become self-shielded 
due to the enhancement in the recombination rate in the IR model. 
Comparing the real and redshift space H I maps we see that in the 
initial stages the ionized and neutral regions respectively appear 
slightly contracted and elongated along the LoS. For the HR model 
the H I distribution exhibits the same behaviour even in the late 
stages. The small H I clumps which appear in the late stages of 
the IR model however, show the opposite behaviour. These neu- 
tral regions appear contracted along the LoS. All the features men- 
tioned above can be understood based on the fact that overdense 
and underdense regions respectively appear contracted and elon- 
gated along the LoS. 

Figure 2 shows the real space power spectrum P r (k) and 
the monopole Po(k) of the redshift space power spectrum for a 
few representative values of xh,. We first discuss the behaviour 
of P r (k) for the HR model. We find that at the very early stage 
of reionization (xh i = 0.9) there is a drop in the amplitude, the 
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shape however remains nearly the same as the dark matter power 
spectrum (xu i = 1.0). At xh i = 0.8, the power spectrum drops 
sharply at k modes smaller than k ~ 0.3 Mpc" 1 , whereas there 
is an increase in power relative to ihi = 0.9 at k modes larger 
than this value. Subsequently, the power spectrum becomes flat at 
low values of xu i (= 0.5, 0.3 and 0.1) and does not retain much 
resemblance to the dark matter power spectrum. In this stage the 
overall amplitude falls with decreasing neutral fraction. Our results 
are consistent with the results reported in the earlier works (e.g. 
Lidz et al. 2007; Choudhury, Haehnelt & Regan 2009; Mesinger et 
al. 201 1). We note that the error-bars in the simulated power spec- 
tra are considerably large at the smallest k value shown in Figure 
2, particularly at low xn j (< 0.3), and it is necessary to be careful 
in interpreting these values. 

The results from the IR model are very similar to those of 
the HR model in the early stage of reionization. Note however that 
small differences can be made out at large k even as early asiHi = 
0.8. These differences are very prominent at the later stages, where 
the predictions of the two models are quite different. The power 
spectrum which is nearly flat for xu i < 0.5 in the HR model, rises 
steeply with slope ~ 1 at k > 0.4 Mpc -1 for xh i = 0.5, 0.3 in 
the IR model. The power spectrum is highly suppressed for xh i = 
0.1, and it shows a power law rise across nearly the entire k range 
that we have considered. Our results are consistent with the findings 
reported earlier by Choudhury, Haehnelt & Regan (2009). 

To understand the evolution of the power spectrum P r (k) 
with xh i in further detail, we study the behaviour of the ratio 
P r /Paa as a function of xh i (Figure 3) at three representative 



k modes k = 0.23, 0.55 and 3.00 Mpc - hitherto referred to as 
large, intermediate and small scales. We could have, in principle, 
chosen a smaller k mode to illustrate the behaviour at large scales. 
The errors however are rather large for k < 0.2 Mpc^ 1 and these 
scales do not provide a very reliable estimate of the behaviour, and 
we need larger simulations to study the behaviour at length scales as 
large as these. We observe that in all the cases the signal falls very 
rapidly as the neutral fraction changes from xu i = 1 to 0.8. This 
decline in the signal is considerably faster than x H ,, which we ex- 
pect from just the decline in the overall H I content. Considering the 
large length scales first, we find that there are three distinct regimes 
in the evolution of P r (k) with xh i. Starting from P r /Paa = 1 
at xh i = 1, it falls to a value very close to 0.1 at xh i = 0.8. The 
ratio subsequently rises to a maxima in the range xh i = 0.4 — 0.5, 
beyond which it again falls for smaller values of ihi- The maxi- 
mum value P r j Paa = 1.7 occurs at xh i = 0.4 in the HR model 
where we have P r / Paa = 0.6 at xh i = 0.1. In contrast we have 
the maximum value P r /Paa — 1.1 at xh i = 0.5, and we have 
P r /Paa f» at Sh i = 0.1 in the IR model. Considering the inter- 
mediate length scales next, for both the HR and IR models the ini- 
tial decline in the ratio does not fall below 0.2. The subsequent peak 
also is relatively much broader compared to the behaviour at large 
scales, and the ratio does not exceed 1 at any value of xh i- We find 
considerable differences between the HR and IR models at small 
length scales. The ratio initially declines to « 0.5 at xh i = 0.9 in 
both the models. In the IR model the ratio subsequently rises to a 
broad maxima (~ 0.7) and then declines to ~ 0.4 at xh i = 0.1. 
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Figure 5. The ratio Pq (fc)/ P r (k) for different neutral fractions and for both the models of reionization. For reference, the horizontal thin dash-dot line (LT) 
shows the value 28/15 predicted by the linear theory of redshift space distortion. 
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Figure 6. Same as Figure 3 but for the ratio Pq / P r . 



The behaviour, however is different in the HR model where the ra- 
tio continues to decline nearly monotonically to w at xh i — 0.1. 

We next consider Pq the monopole component of the redshift 
space power spectrum. The initial (ih i = 1) power spectrum is 
very similar to the real space counterpart, except that the amplitude 
is enhanced. This is approximately consistent with the enhance- 
ment by a factor of 1.87 predicted by linear theory. For both the 
HR and IR models the changes in Pq (k) with xh i are very sim- 
ilar to those seen for the real space counterpart P r (k) discussed 
earlier. We however, notice two exceptions, the first being the be- 
haviour of the power spectrum at ihi = 0.8. The sharp drop at 
k < 0.3 Mpc" 1 seen in P r {k) is not seen in P s (k) which has 
nearly the same shape but lower amplitudes as compared to Pq (A:) 
for iH[ = 1. The other difference is seen at large k (> 3 Mpc -1 ) 
where Po (fc) has a down turn, whereas P r (k) is either flat or an 
increasing function of k. 

Here also we use the ratio Pq /Paa to study the evolution 
with varying xh i, considering the three particular k modes men- 
tioned earlier. The linear theory predicts this ratio to have a value 
1.87 at ih i = 1. We find our results to be consistent with this value 
at the intermediate and large scales shown in Figure 4, the result at 
small scales is slightly larger, possibly due to non-liner effects be- 
coming important at these scales. Considering the large scale first 
for both the HR and IR models we find three distinctly different 
regimes in the behaviour of P S / Paa with varying xh i. In the ini- 
tial stages we find a rapid decline with a minima at xh i = 0.8. 
This is followed by a sharp rise till xh i = 0.4 — 0.5 where the 
behaviour changes and Pq/Paa again declines. The ratio has a 
peak value ~ 1.5 at xh i = 0.4 and falls to 0.6 at xh i = 0.1 



for the HR model, the results for the IR model are some what dif- 
ferent with a peak value of « 0.9 occurring at xh i = 0.5 and it 
falls to « at x H i = 0.1. For both the HR and IR models, the 
xh i dependence shows only two distinct regimes at intermediate 
and small length scales. The ratio falls sharply till xh i = 0.8 and 
subsequently varies smoothly for xh i < 0.8. The ratio P S /Paa 
shows a gradual rise as xh i from 0.8 to 0.6 and subsequently de- 
clines for smaller values of xh i. This rise in Pq /Paa not seen at 
small length scales where the ratio either remains constant or de- 
clines slowly as xh i falls below 0.8. 

We next consider the ratio Pq jP r (Figure 5) which has been 
extensively used in the literature (e.g. Lidz et al. 2007; Mesinger 
et al. 2011; Mao et al. 2012) to quantify the effect of redshift 
space distortion. The earlier works have reported that at large length 
scales this ratio rises to a value greater than 1.87 in the early stages 
of reionization, and then falls to a value slightly less than 1 at later 
stages. Our results show a similar behaviour in both the HR and IR 
models. To study this in more detail we have shown the xh i de- 
pendence of Pq I P r in Figure 6 for the three k modes mentioned 
earlier. We find that at large scales Pq / P r rises to a value ~ 2 
at xh i = 0.8 — 0.9 and then abruptly falls to a value slightly 
less than 1 at xh i = 0.7. The ratio subsequently rises gradually 
and approaches 1 as xh i decreases. The behaviour at intermedi- 
ate scales is somewhat similar, except that we have a sharp peak at 
xh i = 0.9. Also, the ratio again exceeds 1 at xh i < 0.2 for the 
IR model. At small scales we do not find any peak for xh i < 1, 
and the ratio initially declines rapidly to w 1 at xh i = 0.7 and 
declines gradually subsequently. Several earlier works have high- 
lighted the initial rise in Pq jP T at xh i ~ 0.8 as a very prominent 
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Figure 8. Same as Figure 3 but for the ratio P| / Pq . 
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Figure 9. The dimensionless cross-correlation TZ(k) for different neutral fractions and for both the models of reionization. 



feature of the effect of redshift space distortion on the EoR 21 -cm 
signal. We note, that the expected signal itself drops considerably 
at iei = 0.8 and this is possibly not very significant from the 
observational point of view. Further, the sudden rise in Pq jP r is 
possibly because of the rapid decline in P r at ihi = 0.8 which 
appears in the denominator of this ratio. 



Figure 7 shows the ratio P£/Pq which quantifies the 
quadrupole anisotropy of the redshift space H I power spectrum. We 
find that for xh i = 1.0 the results from our simulations are consis- 
tent with the value 1.01 predicted by linear theory at large length 
scales (k < 0.5 Mpc -1 ). At smaller scales the ratio is larger than 
1.01 possibly because of the non-linear redshift space distortion. 
This ratio increases at all length scales for xh i = 0.9. The ratio 
increases further for xh i = 0.8 at large scales (k < 0.3 Mpc -1 ) 
whereas it decreases at smaller length scales. The behaviour of this 
ratio changes significantly atiHi < 0.7 where we find that Pi/Po 
is negative at large scales. The value increases steadily as we move 
to smaller scales, crosses zero at k ~ 1.0 Mpc -1 and is positive 
at smaller length scales. The results for the HR and IR models are 
very similar except that there is a bump at around k ~ 0.5 Mpc -1 
for the IR model. To study the evolution of P% /Pq with xh i in 
more detail we show this in Figure 8 for the three k modes men- 
tioned earlier. Considering the large scale first, the behaviour of the 
HR and IR models are both quite similar. The ratio rises from ~ 1 
at xn i = 1.0 to 1.4 atiHi = 0.8 — 0.9 and then ab-raptly falls 
to m —0.3 at Sh i = 0.7. The ratio remains negative with values 
in the range —0.4 0.6 for smaller values of xh i with the ex- 
ception that it rises to ~ —0.2 for xu i < 0.2 in the IR model. The 
behaviour at intermediate scales is very similar as that at the large 
scales except that the ratio falls to a value in the range 0.2 — 0.4 at 
xh i = 0.7 instead of becoming negative. The ratio declines further 
and is negative (~ —0.1) for xh i < 0.4 in the HR model. In the 
IR model the ratio is negative nowhere and is in the range 0.2 — 0.4 
for xh i < 0.7. At small scales the ratio is constant at ~ 1.2 for 
xh i > 0.9 where after it falls rapidly to 0.8 at xh i = 0.7 and 
subsequently declines gradually to ~ 0.6 and ~ 0.4 in the HR and 
IR models respectively. 



The hexadecapole Pi(k) measured from our simulations has 
very large error bars, and consequently we have not shown these 
here. Larger simulations are required for reliable estimates of 
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Figure 10. This shows A(k) = P X x{k) / P&/±(k) for different neutral 
fractions and for both the models of reionization. The line styles are the 
same as Figure 9. In the left panel the bottom most curve corresponds to 
xh i = 0.9 and the value of A(k) increases monotonically with decreasing 
xh i- 



5 DISCUSSION AND CONCLUSIONS 

We first compare our results with the predictions of the model dis- 
cussed in Section 3 for which P T ', P§ and P£ can be calculated 
using eq.s (9), (10) and (11) respectively. The model requires A(k) 
and IZ(k) which we have directly determined from our simulations 
at different stages of reionization. Several earlier authors (Lidz et 
al, 2007; Mesinger et al, 2011; Mao et al, 2012) have noted that 
this model does not correctly reproduce the real and redshift space 
power spectrum P r and Pq determined from simulations. Our re- 
sults for the xh i dependence of P r /Pa a, Pq/Paa, Po/P r and 
Pi /Pa respectively shown in Figures 3, 4, 6 and 8 confirm that 
the model fails to quantitatively reproduce the results of the sim- 
ulations. We however note that the predictions of this model are 
qualitatively very similar to the results of the simulations, and this 
provides a very useful framework for interpreting our results. In the 
subsequent discussion we have used the framework of this model 
to discuss and qualitatively interpret the xh i dependence of various 
quantities that we have studied earlier. 

Figure 9 and 10 show IZ(k) and A(k) respectively at differ- 
ent stages of reionization. Considering the large scales first we find 
that 1Z ~ —1 at all stages of reionization in both the HR and IR 
models. This indicate that at large scales the distribution of neu- 
tral fraction is anti-correlated with the matter distribution through- 
out reionization, this being a consequence of the fact that the high 
density regions are ionized and only the low density regions re- 
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main neutral. Note that the small neutral clumps produced at the 
later stages in the IR model do not affect this behaviour seen at 
large scales. Assuming that TZ w — 1, the real space power spec- 
trum [eq. (9)] becomes P T / Paa ~ i|,(l - A) 2 . Consider- 
ing the behaviour of A(k) we see that at large scales it starts at 
a value A(k) ~ 0.2 at xh i = 0.9 and gradually increases as 
reionization proceeds, reaching a value ~ 10 at xn i = 0.3. It is 
possible to qualitatively explain the xu i dependence of P r /Paa 
(Figure 3) at large scales using these considerations. We expect 
P T / Paa to decline as A(k) approaches 1 in the early stages of 
reionization, P T / Paa will have a minima as A(k) ~ 1 which 
occurs between xn i = 0.7 and 0.8 in our simulations. The ra- 
tio will subsequently pick up again as A(k) continues to increase 
in the later stages of reionization where P r /PAA ~ x'hi-4 2 - 
This enhancement due to A(k) is compensated by the falling neu- 
tral fraction ihi and P r / Paa reaches a maxima and declines 
as reionization proceeds further. To analyze the redshift space 
power spectrum at large scales using similar considerations we 
have PS /Paa ~ x\ ,[(28/15) - (8/3).4 + A 2 ] and Pi /Paa » 
(4/3)Sh i [(10/7) — A]/ which provide a qualitative understanding 
of the behaviour seen in the simulation. For both of these quantities 
the initial behaviour (A <C 1) is close to the prediction for the dark 
matter distribution. There is a transition in the behaviour around 
A ~ 1 which is followed by a regime (.4 > 1) where the signal is 
dominated by the fluctuations in the H I neutral fraction. We have 
the limiting values (P r , P S , Pi) = »hi-Paa(1, 1.87, 1.90) and 
= i| ,Paa»4 2 (1, 1, -1.3/-4) at A < 1 and A > 1 respectively. 
These, together with the fact that there is a transition between these 
two limiting behaviour at .4 ~ 1 provides a reasonable understand- 
ing of most the results seen at large scales. 

The fluctuations in the neutral fraction A x (k) and the matter 
A(fc) are not perfectly anti-correlated (TZ > —1) at intermediate 
and small scales. The value of TZ increases as we go to smaller 
scales, and the behaviour is similar in the HR and IR models at 
the early stages of reionization (ih i > 0.8). We see differences 
between the HR and IR model at intermediate and small scales at 
the later stages of reionization. The values of TZ are larger in the 
IR model in comparison to the HR model, and TZ becomes positive 
at small scales in the later stages in the IR model whereas it is 
nowhere positive in the HR model. This difference is the outcome 
of the small H I clumps seen in Figure 1 at the later stages in the 
IR model. Considering A(k) also, we find that the behaviour is 
very similar for both the HR and IR models in the early stages, 
whereas significant differences are noticeable in the later stages. 
The increase in A with decreasing ih i is much less pronounced 
at intermediate and small scales in comparison to the large scales 
discussed earlier. It is possible to qualitatively interpret our results 
for the H I power spectrum at intermediate and small scales based 
on the behaviour of TZ and A discussed here. We however do not 
present a detailed discussion of this here. 

It is anticipated that the initial observational attempts (e.g. 
GMRT, LOFAR and MWA) will probe the EoR 21 -cm signal only 
at large scales. The multipole moments of the power spectrum of 
the measured signal hold the key to quantify and interpret it. We 
have considered the monopole and quadrupole moments which re- 
spectively quantify the strength of the signal and the anisotropy 
which arises due to the peculiar velocities. Our simulations show 
that at large scales the signal is expected to fall sharply as xn , de- 
clines from 1 to 0.8 in the initial stages of reionization. The signal 
is subsequently expected to rise to a maxima at ihi ~ 0.5 after 
which it again declines. The prospects of detecting the EoR 21 -cm 
signal are most favourable when the mean neutral fraction is in the 



range xh i = 0.4—0.5. In this range the signal, we see, is character- 
ized by two main features for both the HR and IR models of reion- 
ization. First, the monopole moment of the redshift space H I power 
spectrum Pg(k) is nearly equal to the real space H I power spec- 
trum P r (k), and both of these are comparable to the real space mat- 
ter power spectrum. Second, the quadrupole moment of the redshift 
space H I power spectrum P£(k) is negative with a value which is 
—0.5 times Pg(k). This is in contrast to the value 1.01 predicted 
by linear theory for the ratio P| /Pq of the matter power spectrum. 
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